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Abstract 

We investigate the effect of the Majorana Fermions which are formed at the boundary of a p- 
wave superconductor. When the Majorana overlapping energy is finite we construct the scattering 
matrix S by maping the Majorana zero mode to Fermions for which coherent states are defined and 
a path integral is obtained . The path integral is used to compute the scattering matrix in terms of 
the electrons in the leads . This method is suitable for computing the conductivity. We investigate 
a chiral Majorana Hamiltonian and show that in the absence of vortices the conductivity vanish. 
We compute the conductivity for p wave superconductor coupled to two metallic leads we show 
that when the overlapping energy between the two Majorana fermions is finite the Andreev Crossed 
reflection conductance is finite. 
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1. INTRODUCTION 


At the surface of a topological insulator electrons carry a Berry phase of tt, in the presence of 
an attractive interactions superconductivity is induced. For a low level doping we obtain a p- 
wave topological superconductor. Majorana Fermions appear on the surface of a topological 
insulator in a region where the chemical potential /ie//.(F) changes sign . We consider the 
effect of the Majorana modes on the p — wave superconductor mg. When two metallic 
leads are attached to the superconductor, the Majorana fermion induces resonant Andreev 
reflection or crossed Majorana Andreev reflection. With increasing doping, a regular 
superconductor is expected with the Andreev conductance of the order of (F is the 

tunneling width and A is the superconducting gap), which is much smaller than the Andreev 
conductance carried by the Majorana fermions. The phenomena of Andreev reflection and 
crossed Andreev reflection can be understood from the general properties of the scattering 
S [5Hn]. The scattering matrix S is computed using the continuity equations and the 
unitarity properties. For hnite Majorana energies, it is difficult to obtain the scattering 
matrix S. For such cases it is preferable to represent the scattering matrix S as a Dyson 


series S = T 


.-UAc dt'H,ff.it') 


na, expressed in terms of the leads Hamiltonian. This is 


obtained by integrating the Majorana fermions. This approach has the advantage of taking 
account the renormalization effect for the tunneling matrix element. Experimentally the 
tunneling for the differential conductance is in disagreement with the quantized values [13] . 
The purpose of this paper is to introduce the scattering matrix S as a Dyson 


series S = T 


3 U-CXD dt' Hgff.it') 


Using the scattering matrix we will compute the 
differential conductance for different cases considered in the literatnre The plan 
of the paper is as following: In sec. 2. we formulate the problem in terms of the leads 
and the fermionic representation of the Majorana fermions. In sec. 3. we consider a 
superconducting island deposited on the surface of a three dimensional topological insulator. 
The area outside the superconductor is gaped by a ferromagnetic material. We demonstrate 
that in the absence of vortices the conductivity between the metallic leads vanish. In sec.4. 
we consider two Majorana fermions coupled to two leads and compute the Andreev crossed 
reflection for the p-wave superconductor . Sec .5. is devoted to conclusions. 


2. MAJORANA FERMIONS FOR A P-WAVE SUPERCONDUCTOR 
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We consider a p-wave superconductor described by the ^a{,x,y) the Bogoliubov -de Genes 
fermion operator. At the boundary between the superconductor and the metallic leads 
Majorana zero modes are formed (the chemical potential changes sign). The Bogoliubov -de 
Genes operator contains also the zero modes given by the operator Go(r). The coupling of 
the p-wave superconductor to the two leads is given by, 

Ht = t^ f dy dl{x = -^,y)^^{x = -^,y) + dl{x=^,y)^^{{x=^,y) + h.c. ( 1 ) 

f7=t,4 

dl{x = —^,y),da{x = — f,|/) are the fermions in the left lead and d'l{x = ^,y),d„{x = ^,y) 
represent the fermions in the right lead. 

I the presence of the Majorana fermions we replace the Bogoliubov -de Genes operator 
'^a{x,y) by the zero mode part Co{r^. 

For an even number of Majorana fermions, we replace the zero mode CQ{r) by the repre¬ 
sentation : 

n 

Go(r) = ^2 ^ ^ 72 a-l-F 2 a-l(r) + %aF 2 a{f) ( 2 ) 

a=\ 

The spinors are given by F2a{f) and F2a-i{f) 



(3) 


The two component spinors are localized at the positions r = i? 2 a-i and r = i? 2 a- We 
introduce the fermion operators (I and C,a , a = l,2,3...n. The transformation between the 
two representation is given by: 72 a-i = ^ Cl + Ca > 72 a = ^ Cl - Ca , a = 1, 2, 3...n 
The overlap between different Majorana fermions will introduce the energy for the Ma¬ 
jorana Hamiltonian: H^^^^rana) ^ iea 72 a-i 72 a = ‘^aClCa- This allows to obtain 

the Majorana action. 


= 5^ / * [Cl(*9i)Ca - e.CiC.] (4) 

a=l 


The action in Eq.(4) allows for the construction of the integral for the Majorana fermions 
which will be used for computing the conductivity. For an odd number of Majorama 


Fermions we will have for the 2n -|- 1 Majorana an unpaired Fermionic ,we can choose 
for 7 „+i = W or 7 „+i = ^ Cl+i - Cn+i • 
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3. A CHIRAL MAJORANA FERMION COUPLED TO TWO LEADS 


We consider a grounded superconducting island of radius R deposited on the surface of a 
three dimensional topological insulator. The area outside the superconductor is gaped by 
a ferromagnetic material. We will attach the superconducting island to two leads at 6 = 0 
(left lead) and 9 = n (right lead). We will show that in the absence of vortices the left lead 
is effectively not coupled to the right lead and therefore the conductance vanish. 

3.1. No vortex in the superconductor 

The Hamiltonian at the interface is described by a chiral Majorana Hamiltonian. 

RMaJorana = i (5) 

Jo 

We replace the Majorana fermion 7 ( 0 ) by regular fermions C{6) and C^9) , ^{6) = 
^[C^(9) + C{6)] and expand the fermion in angular momentum states: C{6) = 

The Majorana Hamiltonian takes the Bogoliubov-de Genes form; 

^Majorana = ^ ^ [(C/G - CliC.l)l + {ClC-l - C/Cl;)/ (6) 

l>0 

The Bogoliubov-de Genes eigenvalues for the Hamiltonian in Eq.(6) are: eA=o = 0 and 
eA^o = ^2/. The eigenspinors are ,x(/)( for the zero eigenvalue ) and r]{l) ( for the non zero 
eigenvalues). 

x(0 = ^[c'(i) + c(-/)| 

ri(l) = ^[C(/) = ^[C'(i) - C(-i( 

I 

(7) 


The tunneling Hamiltonian is given by 

Ht = g d\{6 = 0)e*^ — di{6 = 0)e“*^) 7(6* = 0) -I- dl{6 = 7r)e*^ — d2{9 = 7r)e“*^) 7(6 * = tt), 

( 8 ) 


We substitute the eigenspinor ,%(/) and rjJ) and hnd: 
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H, = g{d\{e = 0 )e*^ - d,{e = 0 )e-*^) x{l) + 

i 

g(dl{e = 7r)e*^ - d2{9 = 7r)e"*^)( ^ x (0 “ Y1 

l{even) l{odd) 

(9) 

The Hamiltonian Ht in Eq.(9) is independent from x(/) . The integration of of the Majorana 
fermions in Eqs.(7 — 10) will give a scattering matrix.We hnd that the scattering matrix 
depends only on the right lead! The left lead d\{6 = 0)e*^ — di{9 = 0)e“*^ which 
couples to x(/) will not appear in the scattering matrix! As a result the cross- Andreev 
conductance will vanish. 

3.2. A vortex inside the superconductor 

When a vortex is added to the case considered in case given in 3.1 we need to add the 
impurity Hamiltonian: 

p27T 

Himp. Y d0t{0 

( 10 ) 

^vortexi^ is the Majorana vortex which couple with the strength t( 6 ';ro) to the chiral Ma¬ 
jorana fermion 7 ( 6 *;ho). Due to this coupling t( 6 ';ro) the two leads will be coupled and the 
cross- Andreev conductance will be hnite . The exact result of the Andreev conductance 
will depend on the details of the coupling 7 ( 6 *; ho). 

4. A PAIR OF TWO MAJORANA FERMIONS COUPLED TO TWO LEADS 

We consider a grounded p-wave topological superconductor attached to two leads. Close 
to the leads due to the boundary condition the p-wave superconductor has to Majorana 
modes. We will compute the Crossed Andreev Reflection a process where an incoming 
electron from lead 1 is turned into an outgoing hole in lead 2. In this case a single 
electron at each lead is tunneling into superconductor to form a Cooper pair. 
We consider two half vortices localized in the superconductor at r = i?i = [xi ~ 1 / = 0] 

and r = R 2 = [x 2 ^ = 0]. For this case, we have for the zero modes, 

Co(r) = Cl{r) = 72 [7iFi(f) + 72^2^)] (11) 
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where 71 and 72 are the two Majorana operators. We attach the two leads at [x = — ^, ?/ = 0] 
and [x = ^, 1 / = 0], and dne to the non-locality of the spinors Fi{f) , F 2 (r) the Majorana 
fermions conples to the fermions in the two leads . We consider a sitnation where the two 
Majorana fermions overlap with energy eo- Using the energy eo we constrnct the }j(^o,jorana) 
Hamiltonian. The tnnneling Hamiltonian between the leads and the Majorana fermions is 
given by the Hamiltonian ifp 


Ht = g 

^{Majorana) ^ 


- di{x = -^, 0 ) 


_j7L 

e 2 


7i + 


4(^ = ^,0)e*^2" -d2{x = ^,0)e *‘^2" 


72 


( 12 ) 


d\{x = —^, 0 ), di{x = —f, 0 ) are the fermions in the left lead and dl{x = —^, 0 ), 
d 2 {x = — f, 0 ) are the fermions in the right lead, eo describes the overlapping between 
the two Majorana Fermions (two independent half vortices). The two vortices are localized 
at positions i?i , R 2 and their wave fnnctions are non-orthogonal. We replace the two Majo¬ 


rana Fermions with a single fermion ; 7 i = :^ + C cind 72 = ^ ~ C 

Hamiltonian is given in terms of leads operators and V form: 


The tnnneling 




r'c + cV 


(13) 


The operators and V are expressed in terms of the one dimensional leads: 


ut = 

V = - 


iV2i- 

1 

i\/2 


e^'"idl(x = -^,0) - e ^'*'2 di(x = -^,0) 


+ 


42 L 


e*'? 4 (x = ^, 0 ) - e ^'^ 2 d 2 (x = ^,0) 


;7l 


2dl(x = -^,0)-e ^^2di(x = -^,0) 


^42 


1 r 1 ^ / u „, _i’t 2 , , L 


e* 2 " 4 (x =-, 0 ) - e ^ 2 d 2 {x = -,0) 

(14) 


The action for this case is given by: 


S= dt 


i(C*(*)(* 3 .)C(() + C(*)(«.)(*(*)) - - H,{t) 


(16) 


Using the Grassman integration m (see Eq. (1.191) in Nakahara) for the Majorana Fermions 
4,C we obtain the effective Hamiltonian ife//.(t) for the leads: 

poo 

= (-*4) / dr - r) 

Jo ^ 

V^{t)V{t — T) = (^d[{t)e~^^ — di{t)e^^ + id\{t)e~'‘~^ — id 2 {t)F^'^ ■ 

^ — d\{t — r)e“*^ J- diit — r)e*^ -1- id\{t — r)e“*^ — id2{t — r)e*^ j 

(16) 
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For the electrons in the leads, we nse the right {R) and left (L) movers representation. 
di{x = —f,0) = di{t) and d\{x = —f,0) = d\{t) are the electrons in the left lead (1) and 
d 2 {x = ^,0) = d 2 {t) and d\{x = ^, 0) = d\{t) are the eletrons in the right lead (1). 

d,{t) = (17) 


We apply on the left lead a voltage V and on the right lead a voltage —V/2. As a resnlt, 
we obtain for each lead, two Green’s fnnctions. For the left lead (1) we have Gq^{E,uj) 
(right mover) and Gq^{E,oj) (left mover). 


Gl^^{E,u) 

Gl’^{E,u) 


e{E-f) Q{-E + f) 

u-{E-^)+iO ^ u-{E-^)-iO 

e{-E + f) 0(^-f) 

u + {E - f) + lO^ u + {E - f) - tO) 


(18) 


Similarly, for the right (2) lead we have 


Gl'^{E,u) 

Gl’^{E,u) 


e{E + f) Q{-E-f) 

oj-{E + f)+i0^oj-{E + f)-i0 
Q{-E-f) Q{E+f) 

u + {E+f)+t0^u + {E + f)-t0 


(19) 


Q{x) is the step fnnction which is zero for x < 0 and one for x > 0). The cnrrent in the 
leads is given by: J(x = —^) = en(A"Q’'^ — A^"q’^) = J(x = = ev^N^’^ — Nq'^). v is 

then electron velocity in both leads, is the cnrrent density in the left (1) lead, 

and Nq’^ — Nq’^ is the cnrrent density in the right lead (2) [15]. In order to compnte the 
cnrrent, we will compnte the Green’s fnnctions . The Green’s fnnction will be compnted 
pertnrbatively nsing the effective conpling to the leads Heff.it) given in Eq.(16). Heff.it) is 
represented in terms of one dimensional fermions given in Eq.(17) we have for each lead right 
(i?i, R 2 ) and left (Li, L 2 ) fermions. The pertnrbation theory is controlled by the conpling 
constant . We will compute perturbatively the Green’s function G^’^{E,u; ^), 
G^’^{E,uj] ^) (left leads) and G^’'^(E, cn; —^), G^’^(E, w; —^) (right leads) (the in¬ 
dex 1 and 2 represent the leads and L and R represents the left and right fermions.This 
Green’s fnnction contains the contribntions of the particles-holes, particles-particles, 
and holes-holes in the same and opposite leads. From the Green’s fnnction we extract 
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the self energies for each lead and each mover, and We 

hnd, to order g‘^, the self energies: 


T{uj,Uo) 


—2T(u,Uo)—Ln( 
V V 


-2T(u, Uq)—L n 

V 


t + 

2 

A 

t - 

-^ 

A 

1 + 


A 

1 - 

-2“ ^ 


+ i2T{u, Uq) —sgn(a;) 

V 


i2T(u, cjq) —sgnfo;) 

V 


Fq + (ta + 


( 20 ) 


Where A is the band with, Fq is a damping factor which is induced at high momenta, 
and hwo = eo is the Majorana energy. The imaginary part of the self energy obeys 

and the real part of the self energy obeys = 

gfj^i,R((^) ^ The Green’s funtions are given in terms of the self energies: 


G‘’«(£, c; 5^) = (c - (B - ■' 


G‘yB,a,;T)= + 


il,L| 


U 


-1 


( 21 ) 


The real part of the self energy is used to compute the wave function renormalization function 

Z. 


Z“^ 

f = 




F 

,1 + A 


|a;=0 — Z 
1 


1 _ £F 

2A 


4T(a; = 0,uo)g‘^ 


4g^ 


^'(ro + ^^o) 


( 22 ) 


The tunneling current at the left leads will be given by I{V) = ev{N^’^ — N^’^) (which 
replaces I(y = 0) = en(A"Q’'^ — A"q’'^) the expression for zero voltage) in terms of the 
renormalized Green’s function. 


I(V) = ev(^N^’^ - 




'-A 


r du 

dE / —e 

. 271 


tcjO"*" 


Z0(a;) 


u 


EZ - ^Z 


fFZ 


+ 


Z0(-a;) 


u 


EZ + ^Z + fFZ 



(23) 


Z0(a;) 


u-EZ- ^Z + iTZ 


+ 


ZQ{-oj) 


uj-EZ + ^Z-iVZ^^ 


We will use liV) = ev{^N^'^ — to evaluate the differential conductance for the 

Crossed Andreev reflection Due to the nonlinearity of the effective action, we 

will use the scaling equations [I2l IThl - fTS] for the coupling constant g^. The scaling of 
determines the width T. We And the Renormalization Group equation for the width T, 
^ = —const.T"^ with I = log . The solution r(R) as a function of r(R = 0) is given by: 


f(R) = 


f(R = 0) 


(1 + const.Log[^]y 


(24) 


This solution will be used in Eq.(23), where f is replaced by f (R). Substituting f (R) gives 
us the result for the differential conductance 


dI{V) 

dV 


e 

h 


-AZ 




'-AZ 


TZ 


f/FD (D + —Z) + (D - —Z' 
27r L((]_e)2 + (rZ)2JdDV^'^-^-^ 2 ^ 2 ' 


eR. 


^ ^ / de^ 


h 


'-A 


f(R) 


+ 


f(R) 


27rLr^_e)2 + (r(R))2 (^ + e)2 + (r(R))2J 






(25) 


We And that for a pair of vortices the Andreev crossed reAectioion obeys 


IV-5-O 


h 


Figure (1) shows shows the diAerential jy vconductivity for the Andreev crossed reAection 
as a function of the voltage diAerence between the two leads. We observe that in the limit 
R —> 0 the maximum value for the conductance is obtained. This result follows from the 
scaling equation for the width given in Eq.(24). 

Comparing the diferential conductivity with the experiments [13] one observes that the 
perfect quantization is not achieved this suggest the possibility that the width is controlled 
by additional operators causing r(R) not to Aow to zero when R —)■ 0. 


5. CONCLUSIONS 


In this paper we have introduced a new method for computing the conductance in the 
presence of the Majorana Fermions. We mapp the problem of Majorana Fermions to regular 
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FIG. 1: jy the differential conductivity for the Andreevv crossed reflection 


Fermions for which a path integral and the Berry phase are obtained. This allows us to 
integrate out the Majorana Fermions and allows us to obtain the scattering matrix S as a 


Dyson series S = T 


3 D-oo ^ef 


. Using this method we have compute the differential 


conductance for different cases ,Achiral Majorana Fermion coupled to leads with and without 
an additional vortex and studied the Andreev crossed reflection for a pair of Majorana 
coupled to two leads. We have computed the differential vconductivity for the Andreev 
crossed reflection as a function of the voltage difference between the two leads. We observe 
that in the limit U —)■ 0 the conductance reaches the maximum value. 
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